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Abstract 

We investigate the classical integrability of the Alday-Arutyunov-Frolov model, and show 
that the Lax connection can be reduced to a simpler 2x2 representation. Based on this result, we 
calculate the algebra between the L-operators and find that it has a highly non-ultralocal form. 
We then employ and make a suitable generalization of the regularization technique proposed 
by Maillet for a simpler class of non-ultralocal models, and find the corresponding r- and s- 
matrices. We also make a connection between the operator-regularization method proposed 
earlier for the quantum case, and the Maillet's symmetric limit regularization prescription used 
for non-ultralocal algebras in the classical theory. 
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1 Introduction 

There has been a renewed interest in the last decade to the theory of integrable models and 
its various techniques, due to their numerous applications in the understanding of the AdS/CFT 
correspondence (for a comprehensive review see [1,2]). On the other hand, this study has stimulated 
a deeper investigation of the subtleties associated to the quantization of integrable systems. The 
standard methods which work for the simpler classical integrable models [3-6] often fail when 
considering a more complex theory such as the string theory on AdS$ x S" 5 background, or even its 
smaller subsectors . 

In this paper we address one such subsector, the su(l, 1) subsector, and point out some of the 
new interesting features and problems that arise in the resulting Alday-Arutyunov-Frolov [AAF) 
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model [7-13]. The first attempt to probe its quantum integrability has been done in [14], where it 
has been shown that the quantum integrability of the AAF model can indeed be verified directly 
by checking the S'-matrix factorization property, which is a necessary condition for the quantum 
integrability of the system. However, the technically complex perturbative calculation was possible 
to carry out only up to the 1-loop order, and it does not seem to be possible to generalize the 
perturbative test of this property to all-loop order. Besides, as we have emphasized in [14], the 
perturbative calculations for the AAF model contain several subtleties and fine points. Moreover, 
the most interesting properties of the integrable systems are generally non-perturbative in nature, 
and, therefore, developing a more strict approach is very desirable. 

Thus, our goal in this paper is to develop the more reliable inverse scattering method for the AAF 
model, which does not utilize any perturbative calculations. The standard path to quantizing an 
integrable system is to start from the classical integrability and the corresponding Lax connection, 
and find the algebra of the corresponding monodromy matrices, which determines the classical 
r-matrix of the model. The quantization is then usually achieved by putting the system onto a 
lattice, with the use of the classical algebra. While this works for a number of classical models, it 
does not work in such a straightforward manner for many interesting models, due to the following 
several reasons. First, while the classical algebra may be simple, finding the lattice version of 
the Lax pair is a very non-trivial problem. Even for the simplest classical models, for example, 
the sine-Gordon model, the general techniques of constructing the lattice counterparts do not give 
a simple local Hamiltonian, even though they may exist in principle [15-18]. In some cases, for 
example, the non-linear Schrodinger (NLS) model, one may write the corresponding quantum 
equations directly in the continuous case without any evident problems. This, however, fails for 
other models, as we have shown in [19,20] for the case of the Landau-Lifshitz (LL) model [21], for 
which the direct generalization of the classical equations to the quantum ones leads to meaningless 
singular expressions, due to the ill-defined operator product in the same point. As was shown in 
[19,20], a more mathematically correct procedure to obtain the corresponding non-singular quantum 
expressions is to introduce a special operator regularization, and necessarily construct the self- 
adjoint extensions. Together with the self-adjointness of the quantum Hamiltonian this was shown 
to reproduce both the correct spectrum and the factorization property of the S'-matrix. The root 
of the problem in the LL model was the very singular S"(x) type of interaction in the quantum- 
mechanical description. We emphasize that even in the NLS model one in principle should construct 
the corresponding self-adjoint extensions, although, in this case, the result is non-essential for the 
integrability of the system. 

The AAF model, which contains only fermionic degrees of freedom, is another interesting ex- 
ample of such singular theories. Here one also has to deal with the 5"{x) type of potential, and, 
therefore, all the methods and operator regularization techniques developed for the LL model could 
also be applied in this case. However, unlike the LL model, there is a further complication in the 
AAF model, apparent already on the classical level. As we will show below, the algebra of the L- 
operators has a non-ultralocal form, which essentially prevents one from using the standard methods 
of the integrable systems and quantize the model via the Bethe Ansatz techniques. Although the 
non-ultralocality appears in many interesting models, e.g. principal chiral mode, 2d gravity etc., 
the non-ultralocality in the algebra of the .L-operators is usually exhibited in terms proportional 
to d x 5{x — y). However, the corresponding algebra of the L-operators for the AAF-model has a 
surprisingly more complicated structure, containing higher order non-ultralocal terms of the type 

Thus, the AAF model is an interesting model in this context, which exhibits both difficulties 
simultaneously - the singularity of the interaction, and the non-ultralocality of the algebra of the 
L operators. The singular nature can be dealt with in the same manner as it was done for the 
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LL model [19, 20]. However, handling the non-ultralocal nature of the algebra is not an easy 
task. There do not exist any satisfactory nor standard methods to deal with such algebras. The 
main prescriptions in this direction are due to Maillet [22-25] and, alternatively, to Faddeev and 
Reshetikhin [26]. The latter method, being more elegant and physically clear, is, nevertheless, hard 
to use in practice for more involved models, and, moreover, it still requires putting the system on 
the lattice. 1 The method due to Maillet, however, does not use any lattice regularization, and 
although it is not obvious how to quantize such systems, some essential progress in understanding 
the integrability of such models, e.g., the complex sine-Gordon model and non-linear sigma models, 
has been made in the classical theory (for more recent applications see [29-33] and the references 
therein). 

One of the main issues regarding the method proposed by Maillet is the ad hoc construction 
of the symmetric limit procedure, in order to obtain well-defined algebras between the monodromy 
matrices. 2 In this paper we will argue that such symmetric limit procedure is the result of the 
regularized operator product in the quantum theory, and which naturally appears when one takes 
the classical limit H — >■ 0. 

Another interesting result we have obtained is that the AAF model admits a 2 x 2 Lax pair 
representation. Let us remind, that in [7] it was shown that one can obtain a 4 X 4 representation, 
starting from the Lax pair for the full superstring on AdS§ x 5 5 and carefully eliminating the 
degrees of freedom, in order to obtain the su(l, 1) subsector. Moreover, it was proved there that 
the the zero-curvature condition is satisfied on the equations of motion. Here we demonstrate also 
the inverse, namely, we show that one can derive the equations of motions from the zero-curvature 
condition. In the process, by analysing all the constraints and the resulting equations, we show 
that the Lax connections can be recast in the 2x2 matrix form, which essentially simplifies the 
computation of the algebra of the L-operators. This is quite remarkable, since the much simpler 
fermionic Thirring model admits only a 3 x 3 matrix representation for the Lax operators [35,36], 
which makes it complicated to use the Bethe Ansatz to find the spectrum. Surprisingly, the AAF 
model appears to be simpler in this sense. 

We also give the complete account of the non-ultralocal algebra between the L-operators, which 
is hard to compute explicitly due to a very complex Dirac bracket structure between the fermionic 
fields. We note here, that as has been shown in [14], in the process of checking the S-matrix fac- 
torization property via perturbative 1-loop calculations, we have discovered some missed numerical 
factors in the AAF Lagrangian in the earlier works [7, 9] , which essentially changed some of the 
results. Here, we also correct the missed factors in the Diract bracket structure of [7]. We show 
that in order to describe the integrable structure of the AAF model, one needs to introduce three 
independent matrices r, s\ and S2- The latter is due to the additional higher order non-ultralocal 
term in the algebra. Moreover, we derive the algebra for the transition matrices and show that it 
has exactly the same form as the Maillet algebra with the non-ultralocality of the simpler type, 
containing only the terms proportional to d x 5(x — y). Thus, we show that the effect of the higher 
order non-ultralocality can be absorbed into only two independent matrices, the modified r and s 
pair. 

Our paper is organized as follows: In Section 2, we set up our notations and analyse the 4x4 
matrix representation in details. We show that the AAF equations of motion follow from the off- 
diagonal part of the zero-curvature condition, while the diagonal part gives some highly non-trivial 
constraints, satisfied on the equations of motion. In Section 3, we make the crucial observation that 

For a recent attempt to apply this method to the strings on AdSs x S see [27, 28] 

2 Remarkably, in some cases, such as the 2d gravity coupled to a dilaton field, the algebra is well-defined (in the 
infinite space limit), despite the non-ultralocality in the algebra of the L-operators, due to the presence of the dilaton 
field and its assymptotic behavior [34]. 
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only half of the set of all equations, that follow from the zero-curvature condition, are indeed inde- 
pendent, and show that due to this doubling of the equations, one can reduce the Lax connections 
to a 2 x 2 matrix form. In Section 4, we briefly explain the Faddeev-Jackiw procedure to find the 
Dirac brackets, and give the corrected canonical structure between the fermionic fields. In Section 5, 
we explicitly calculate the algebra between the L operators and show its non-ultralocal nature. As 
an interesting consequence, we find that the field-independent truncation of the algebra corresponds 
to the fermionic version of the Wadati model, which can serve as a more simple characteristic ex- 
ample to analyse non-ultralocal algebras. In Section 6, we recapitulate the Maillet symmetric limit 
procedure to deal with such algebras, and give its generalization adapted to the more general case 
of an algebra containing also terms proportional to d^.5(x — y). We also briefly explain the oper- 
ator product regularization method proposed in [19,20], and argue that Maillet 's symmetric limit 
procedure appears naturally in the classical limit of the regularized quantum case. In conclusion 
we discuss some open problems and the future work. Finally, we collect some important technical 
details in the appendices. 

2 Alday-Arutyunov-Frolov model: 4x4 Lax connection 

In this section we setup our notations and give the complete analysis of the classical integrability, 
using the 4x4 Lax representation originally found in [7]. Let us note, that it was claimed there 
that the corresponding zero-curvature condition is satisfied upon the substitution of the equations 
of motion. Here we also prove the inverse: the equations of motion of the AAF model follow from 
the zero-curvature condition in the 4x4 representation. This is a very non-trivial result, which we 
present in detail below. In the process we will show that the resulting equations and constraints 
are such that one can reduce the Lax representation to a 2 x 2 form. 

The AAF model is obtained (for a complete analysis, see the original paper [7]) by starting from 
the full superstring theory on AdS^ x S 5 and consistently reducing it to the su(l, 1) subsector. The 
remarkable characteristic feature of this truncation is the elimination of all the bosonic degrees of 
freedom through the constraints. Our starting point is the Lagrangian of the AAF model in the 
form (see appendix A for our notations): 



Here, as explained in [14], we have introduced two coupling constants g<± and g%. The main result of 
the analysis in [14] was the necessary relation between the coupling constants 52 and g% in order the 
guarantee the quantum integrability of the model. Namely, in order for the S'-matrix factorization 
property to hold, up to the 1-loop order, the following relation must be satisfied: 



We will show below that the same condition must also hold for the classically integrable theory. 
For now, we will consider a more general theory defined by (2.1) where the constants 52 arid 53 are 
independent, and below we will show, that the constraint (2.2) should be imposed already in the 
classical theory, from the condition of classical integrability. 

To analyse the classical integrability, it is convenient to write the equations of motion, following 
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from the Lagrangian (2.1), for each component ipi and ip2 separately: 



o 



(2.3) 



(2.4) 



This form, however, is still not suitable to proceed with the analysis, since the time derivative of 
the spinor components enters also into the cubic and higher order terms in (2.3) and (2.4). To this 
end, one can substitute in these higher order terms the time derivatives of the spinor components 
upon multiple usage of the equations (2.3) and (2.4), until the higher order terms depend only on 
space derivatives, and the equations of motion have the form: 



We stress that the convergence of this procedure is guaranteed by the fermionic nature of the fields. 
After very lengthy and tedious calculations, one obtains the expression (B.l) and (B.2) of appendix 
B. Even though the resulting equations are rather cumbersome and have a more complicated form 
in comparison to the original equations (2.3) and (2.4), containing terms up to the seventh order 
in the fermions and their space derivative, they will greatly simplify obtaining a number of very 
non-trivial relations, which will play a central role in establishing the classical integrability. 

Let us now turn to the Lax representation. The 4x4 representation was given in [7], and there 
it was shown that the zero-curvature condition is satisfied upon the substitution of the equations of 
motion. In order to establish the inverse, namely, that the zero-curvature condition produces the 
equations of motion (2.3) and (2.4), we will first generalize the construction of [7], write down the 
equations following from the zero-curvature condition, and find the general conditions upon which 
the off-diagonal terms produce the equations of motion, while the diagonal terms give some very 
non-trivial identities. 

The 4x4 Lax connection has the form: 3 



where Iq = cr 3 <8> cr 3 and Jo = c 3 <8> 1, and the other quantities in (2.6) and (2.7) are defined as 

3 As explained in [7], this Lax representation was obtained by starting from the one for the full superstring on 
AdS$ x S 5 and consistently reducing it to the su(l, 1) subsector. We also note, that under the scale transformation: 
ipi — > api, the Lagrangian (2.1) transforms as: Jf(g2,gs) — > c 2 J? (g' 2 , g's), where g' 2 = <?gi and 513 = c 4 p3. One can 
use such transformation to set: g 2 = 1. The Lax connection (2.6) and (2.7) is written after making such a scaling. 




(2.5) 



L = e r) i +^ T) Jo+^ 

Li =^ <T) /o + ^ <T) ^o + A, 



(2.6) 
(2.7) 
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follows: 1 



At) 
SO 



-(1 + ^)(^-^) + |VpV, 



CI 



i 



(i + V^W'-V'V), 



4VA 



2J 



(2.8) 

(2.9) 
(2.10) 

(2.11) 

(2.12) 
(2.13) 

where the functions Zj = ZjQu), i = 0, . . . ,4, given explicitly in appendix A, depend on the spectral 
parameter fi, and: 



Furthermore, the off-diagonal matrices A T and A CT have the following form: 

At = [ir Jo, ho + ij] - d T (i 3 e - he) , 
a cr = k Jo, z 3 + ue] - d a (i 3 e - he) , 



i 



7r 



V'/o -;/' - V^V + 2^^ - 4^ 



(2.14) 
(2.15) 



Finally, the matrices and have the form: 



a 



( j/i \ 

v 2 

i/ 3 

\ i/ 4 o o o / 



/ 











^1 


\ 










z/ 2 












'h 










V 


v\ 














with OQ = l + ~(i/nl>). (2.16) 



The fermionic degrees u m ,u m , m = 1,...,4 will be determined below from the zero-curvature 
condition: 

d Q L l -d 1 L Q -[L Q ,L 1 ]=Q. (2.17) 

In order to write the complete set of equations following from (2.17) in a more compact form, 
it is convenient to write the A T and A CT matrices in the basis a % (g) a 3 , where a 1 are the usual Pauli 
matrices. It is easy to check the following relations: 



where 



e = aoa 1 ® cr 5 'xij and e = ia$o % ® a' J 'xij , 



i (Vi + V2 + V3 + 1/4) , X12 = I (^1 - Z/2 + ^3 - ^4) , 

i 1 



(2.18) 



X21 

X3i = Xi3 = 0. 



(2.19) 



We write these quantities explicitly in terms of the spinor components ipi and ip2 in appendix G. 
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Similarly, one obtains the relations between Xij anci Finally, introducing: 

Vif = hXij ±ihXij, Aii" ) = -/4i' ) and (2-20) 
one arrives at the desired form: 

A T = A;V®^' with A$ = 2ia 0lr ^ ] - d T (cw4 7) ) > ( 2 - 21 ) 
A„ = A^V ® cr 3 ' with At J ' = 2ia 7 CT 4 +) - 5 CT (a ^7 } ) • ( 2 - 22 ) 

Writing the zero-curvature condition for the diagonal and off-diagonal parts separately, we find the 
following equations: 

do^Io + d a{ a) Jo - d^Io - d^Jo - [A T , A(j] = 0, (2.23) 
aoA, - (?i A T - eS r) [Jo, K] - A a) [Jo, K] = 0. (2.24) 

We first analyse the equation (2.24) for the off-diagonal part of the zero-curvature condition 
(2.17), and show that the equations of motion (2.3), (2.4) follow from it. Using the expressions 
(2.21) and (2.22) for the tensors At and A^ , one can write the expression (2.24) in the component 
form: 

doA'J - diA!? - 2ie 3ki A k JHi ] + 2ie zkl K^^ ] = 0, for i, j = 1, 2. (2.25) 

Substituting all possible values for the indices and using the explicit expressions (2.8) - (2.11), 

one obtains the following system of equations: 

hM ( i j) + UM<£ 3) + hhM^i ) + hhM[f + hhM^ + hUM { ^ = 0. (2.26) 

The explicit form of the functions Mg , M4 , P , and M^f 1 are given in appendix 

C. By using the explicit dependence of the li functions on the spectral parameter fj,, given in the 
equation (A. 3), and analysing its various values and asymptotics, one can show that the equation 
(2.26) is equivalent to the following set of constraints: 

(2.27) 
(2.28) 
(2.29) 
(2.30) 

Despite the complicated dependence of the M^' functions on the fields, one can show that by an 
appropriate choice of V{ in (2.16), the above system of equations reproduces the equations of motion 
(2.3) and (2.4). This in turn will imply that the equations (2.27) - (2.30) are not independent, 
leading to the construction of the 2x2 Lax connection, which we present in the next section. 

We start with the equation (2.27). Using the formulas in appendix C, it is easy to show that, 
for any choice of the indices, the equation (2.27) is equivalent to the following constraint: 

3i7r " d j a = 0, (2.31) 

or, more explicitly: 

V4>i + V4>2 - fati. - r 2 4>2 = -idi {W) , (2.32) 



M 3 (y) 


+ M[? 


= 0, 




-mm 


= 0, 


mW 


- mW 


= 0, 




- M[f 


= 0. 
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where we have used the 7 T and 7 CT defined correspondingly in (2.15) and (2.14). The dependence of 
the fields Ui, on the fermionic fields must be chosen so that this constraint is satisfied. 

Let us now turn to the equation (2.29), from which the fields Ui, V{ are determined. Considering 
all possible choices of indices using the formulas in appendix C, as well as inverting the 

relations (2.19), one arrives at the following equations: 

+ -^(aojrh + ^-^(d (a D k ) = 0, for k = 1,2, (2.33) 

and 



A 2VA 



2i^ T di(aov m ) - 2i'y a d (a i' rn ) + -^^7 CT a 7 CT i>m - ^-j^7adi(a i> m ) 

CaoTr^mH 7^C,d G {a u m ) = 0, for m = 3,4, (2.34) 



A 2VA 
where have we denoted: 

C := 2J - ^ (^V - ^pV) " (2.35) 

Thus, we see that the equations following from the off-diagonal part of the zero-curvature condition 
are not independent, and there are in fact only two independent equations. 

Before proceeding to the analysis of the equations (2.33) and (2.34) in all orders in ipi, one can 
first readily compare the linear terms with the ones arising from the equations of motion of the 
AAF model (2.3) and (2.4). One of the possible choices is the following:' 3 

One may then examine the connections between the equations (2.27) - (2.30). The explicit relations 
are given in details in appendix C, and we conclude that there are only two independent equa- 
tions, following from the off-diagonal part of the zero-curvature condition, which we take to be the 
following: 

M^ j) + M^p = 0, (2.37) 
M W ~ M W = 0. (2.38) 

Substituting the formulas (2.8) - (2.15) into (2.38), and considering now all orders in the fields, one 
can derive the dynamical equations (B.3) and (B.4) for ip\ and i[>2 following from the off-diagonal 
part of the zero-curvature condition. 

One immediately sees that these equations do not seem to coincide with the equations of motion 
for the AAF model (2.3) and (2.4). This is, however, due to the presence of the time derivatives of 
the spinor components in the higher order terms. As we discussed earlier, one can, upon multiple 
substitutions of the time derivatives of the spinor components into the cubic and the higher order 
terms, eliminate such dependences, and arrive at the equations of motions where all the cubic and 

5 We also note, that this choice is consistent with the involution, under which ~ Ci 1* ~ ~7r, 7^ = — 7o- and 
do = a . 
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higher order terms depend only on the fields and their space derivatives. After this very lengthy 
and tedious elimination procedure one arrives exactly at the equations (B.l) and (B.2), provided 
the constraint (2.2) is satisfied. 6 Now, using these equations, it is simple to check the remaining 
independent equation of the off-diagonal part (2.37), which, as we discussed above, is equivalent to 
the constraint (2.32). 

Let us now consider the equations arising from the diagonal part (2.23) of the zero-curvature 
condition (2.17). Writing the A T and A CT matrices in the form: 



( 
















A? 








A? 








\K 








o > 



A T 



one easily obtains the following constraints: 



and A n 












Ai \ 








A 2 








A* 








K 








0. J 



(2.39) 



0, 
1 
2 



+ 22 = 0, 



where we have denoted 



[K,K] = diag(<t> u ,<fr",tp u ,<T)- 



(2.40) 
(2.41) 
(2.42) 

(2.43) 



It is important to stress that, as it was the case when considering the off-diagonal part of the zero- 
curvature condition, each of these equations appears twice. Thus, we come to the conclusion that 
each independent equation that follows from the 4x4 representation of the Lax connection appears 
exactly twice in the full set of equations and constraints. This important observation will lead us 
to the construction of the 2x2 Lax connection in the next section. 

The verification of the equations (2.40), (2.41) and (2.42) again requires very lengthy calcula- 
tions, and the usage of the equations (B.l) and (B.2). 7 Nevertheless, it is quite remarkable that 
these equations are indeed satisfied. In appendix D we give some useful expressions and additional 
technical details related to the above constraints. 

Let us also address one subtlety which we have so far ignored. Since the 4x4 representation 
for the Lax connection was obtained from the original Lax connection for the full superstring on 
AdSs x S 5 , and, as it is well known, there is no matrix representation for the psu(2, 2|4) superalgebra, 
there is a possible anomaly in the diagonal part (for a detailed discussion, see [12]). In other words, 
the diagonal part of the zero-curvature condition for the full superstring on AdS§ x S* 5 should be 
generalized to be equal to a term of the form T4 X 4(^)14, where T^^ijji) is some function depending 
on the fields and I4 is the 4x4 unit matrix. It does not vanish, in general, when reduced to smaller 
subsectors, which is, for example, the case for the reduction to the AdS% x S 3 subsector [12]. For the 
AAF model, however, the equation (2.42) is essentially the condition that the anomaly T^x^ipi) 
vanishes. Indeed, it is easy to see that the anomalous term T^xA^i) can be written as follows: 



+ 



22 



) 



(2.44) 



6 In terms of the rescaled fields (see footnote 3), this constraint becomes: g\ = g-j, = 1. 

7 We emphasize that the equations (B.l) and (B.2) have already been obtained from the off-diagonal part of the 
zero-curvature condition (2.17), and, therefore, we are allowed to use them. 
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It is interesting to note, that for the 2x2 representation of the Lax connection, which we will give in 
the next section, the anomaly T2x2(ifti) identically coincides with the anomaly T^^ipi). Therefore, 
the anomalous term vanishes in both cases. 

Finally, it is tempting to try to generalize the construction of the Lax connection (2.6) - (2.16) 
in such a way that the classical integrability holds without imposing the quantum constraint (2.2), 
which was originally obtained in [14] from the 5-matrix factorization property in the 1-loop order. 
Without giving here the explicit analysis, which is quite tedious and follows the same type of steps 
we outlined above, we can state that under no deformation of the parameters, or addition of higher 
order terms in the formulas (2.8) - (2.16) the AAF model is classically integrable for arbitrary 
constants g 2 and g^. Hence, the quantum constraint (2.2) must also be imposed in the classical 
theory. 

Thus, we have strictly proved the classical integrability of the AAF model, based on the 4x4 
Lax connection, provided the constraint (2.2). 



3 2x2 Lax connection 



As we showed in the previous section, the equations that follow from the zero-curvature condition 
(2.17) and the 4x4 Lax connection are not independent, and each independent equation appears 
exactly twice in the set of all equations. This simple observation allows one to reduce the Lax 
connection to a 2 x 2 representation. Indeed, the 2x2 Lax connection can be written in the 
following form: 



£o = d r) /o + £i r) Jo + A T , 
A=^ CT) /o + rf CT) Jo + A CT . 



(3.1) 
(3.2) 



Here ^\ £p°^ and £^ are defined by the same formulas (2.8) - (2.11), and the 2x2 matrices 
Io and Jo have now the form: 

Io = m^2 and Jo = rj 2 cr 3 , (3.3) 

where r]i and r\2 are some constants which we will fix below. The off-diagonal matrices A T and A CT 
now take the following form: 



A r 

Acr 



d T [he - ue) , 
dJhe-iJ), 



where the matrices 9 and 9 are defined by: 



ao 



ui 
v 2 



and 9 = iao 



h 



(3.4) 
(3.5) 

(3.6) 



As before, we will determine the fermionic degrees v m and u m , m = 1, 2 from the zero-curvature 
condition: 

d £ 1 -di£ o -[£ ,£ 1 ]=0. (3.7) 

Introducing the notations: 



Cl = \ ("1 + V2) , C2 



(ui - u 2 ) 

(pi - h) 



(3.8) 
(3.9) 
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from which we construct the combinations: 



P^=kCi±iUCi, M +) = "4 +) and ~4 +) =A + \ (3-10) 



and writing: 



A T = AX with A T = 2i m a ~fM +) - d T (a /4 } ) , (3.11) 
A (T = A> i with A; = 2ir ?2 a 7^! +) - a CT (a ^" } ) , (3.12) 

it is easy to show that the zero-curvature condition (3.7) can be written as the equations for the 
diagonal and off-diagonal parts, similar to (2.23) and (2.24): 

(d<4 CT) - d4 r) ) Jo + (d^ ] - di£i T) ) Jo - [A T , Aa] = 0, (3.13) 

d A a - dtK - Ci r) [J 0l A,] - [Jo, K] = 0. (3.14) 

One can then analyse these equations similarly to the 4x4 representation case. For the off- 
diagonal equation (3.14) one can show, repeating each step of the calculation in the previous section, 
that the resulting equations coincide with the equations (2.37) and (2. 38), 8 provided the constraint 
(2.2) on the coupling constants <? 2 and g%, as well as the relations: 

r/2 = 1, v\ = ip2, V2 = V>2, V\ = —ipi and z? 2 = tp*. (3.15) 
To simplify the analysis of the diagonal part (3.13), we write: 

A T =( » A „') and A„=( » ^). (3.16) 

Noting that the second term in (3.13) is equal to zero, due to the identity (2.40), and denoting: 

[A r ,A CT ] =diag(<t>,<f>), (3.17) 
it is easy to see that the off-diagonal part (3.13) reduces to the following equation: 

m (d & ] - d4 r) ) -<p = r 2x2 (vO, (3.18) 

where r 2 x2(V'j) is the term that may appear due to the possible anomaly in the diagonal part, as 
we discussed in the end of the previous section. Using the relations (3.15), which were determined 
from the off-diagonal part, and the identity (2.41), one can show, that: 

<t> = 4> {11) , (3.19) 



(V>*) = Ihm-q^-m^]- ( 3 - 2 °) 



Remarkably, it follows from (2.44) and (3.20) that the anomaly r 2x2 (V'i) m the 2x2 case coincides 
identically with the anomaly F^x^V'i) of the 4x4 case, provided 771 = 1. Thus, in both cases the 
anomalous terms vanish. 

8 We stress that the equations of motion for the AAF model (2.3) and (2.4) follow from the equations (2.37) and 
(2.38). 
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4 Dirac brackets 



Having derived the reduced 2x2 Lax representation for the AAF model, we turn our attention to 
its Poisson structure. This analysis has already been carried out in the original paper [7], however, 
due to some missed coefficients in their original Lagrangian, we carefully redo this derivation. Our 
first step is to obtain a Hamiltonian formulation of the model. It is convenient to rescale the fields 
in the Lagrangian (2.1) as follows: if) — > J~2ip. Then, the rescaled Lagrangian becomes: 9 

jsf = - j - 1 (ipp d ip - doV^V) + i~ (V^ 1 ^ - + W> 

+ ^7? ^ ^p" 9 ^ ~ W 9 $p^) - ^Qji^ i^f d ^p 5 d^. (4.1) 

The AAF Hamiltonian can be obtained by the standard Legendre transform: 



where the canonical conjugate momenta are defined by: 



dtp 



<9«£? i V\g<2 , 5 5 V\g3 

= — = -p> + -jj^- {-if) d!tl)p 1p + d X lf)lf) p°ip) - ^J3" 



so that the Hamiltonian becomes: 



(4.2) 

-^(^)V^, (4.4) 
ipif). (4.5) 



Since the Lagrangian (4.1) is linear in the time derivatives, the canonical momenta are indepen- 
dent of the time derivatives of the fields. Hence, the attempt to eliminate the time derivatives of the 
fields in favour of the canonical momenta fails, and one must analyse the constrains in the theory 
and construct the Dirac brackets [37]. This is rather a cumbersome procedure for the AAF model, 
which can be avoided by utilizing the equivalent prescription by Faddeev and Jackiw [38,39], which 
we briefly review below. 

4.1 Overview of the Faddeev- Jackiw formalism 

Faddeev and Jackiw's method is based on the observation that a conventional second order in 
time derivatives Lagrangian can always be converted to a first order in time derivatives Lagrangian 
by the exact same Legendre transform used to go from the Lagrangian to the Hamiltonian formu- 
lation. Let us briefly describe this approach by considering a general first order in time derivatives 
Lagrangian: 10 

L = ai (Oti-V(0, (4.6) 

9 It is interesting to keep the coupling constants gi and gs independent throughout the analysis carried out in this 
section and, therefore, uncover the explicit dependence of the Dirac brackets on them. However, when considering 
the algebra of Lax operators in section 5, we shall use the rescaled fields in which the constraint is g\ = gz = 1. 

10 The Lagrangian considered in this example describes a typical mechanical system. The generalization to a field 
theoretical setting with the inclusion of anticommuting variables is straightforward. It will be considered in detail for 
the AAF model in the next section. 



13 



where £j denote the 2n phase space coordinates: 



£i=Pi, i = l,...,n and & = qi , i = n + 1, . . . , 2n, 

with the sum over repeated indices, as usual, implied, and where aj(£) is an arbitrary function of 
the £j, without explicit time dependence. Noting the absence of first order time derivatives in the 
combination: 

- L, 

when defining a Hamiltonian by the standard Legendre transform, it is possible to identify the 
potential V with the Hamiltonian: 

dL 

H = - f £ i -L = V. (4.7) 

Thus the first term on the right hand side of (4.6) defines the canonical 1-form: a(£) = 
The Euler-Lagrange equations obtained from (4.6) have the form: 

= wth WB = _^__^. (4.8) 

If the 2-form: uj = da = d£i d£j is nonsingular, then the matrix Wjj is invertible, and (4.8) can 
be recast in the form: 

dV 

(4 - 9) 

Since V is the Hamiltonian for Lagrangian (4.6), the equations (4.9) are also Hamiltonian: 

dV 

provided one defines the bracket such that: 

It is important to emphasize that it was not necessary to consider any constraints in the analysis 
so far. They only appear in the case where the matrix cuij is singular, and a more involved analysis 
is required. This discussion is, however, out of the scope of the present work, since, in the case of 
the AAF model, the matrix is invertible. We refer the interested reader to [38,39]. Finally, we 
stress that the bracket (4.11) coincides with the one obtained through the Dirac procedure [40]. 

4.2 Faddeev-Jackiw formalism for the AAF model 

In this section we apply the prescription due to Faddeev and Jackiw to the AAF model. We 
start by noting that the Lagrangian (4.1) admits the following decomposition: 

if = if fcin -^, (4.12) 

where Jf? is the AAF Hamiltonian written in equation (4.5) and J^kin stands for the kinetic part of 
the Lagrangian: 

-^-e^{^fd a ^p 5 d^. (4.13) 
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In order to obtain an analogue of (4.6), we must extract the canonical 1-form from =£?fcj n , namely, 
we must write: 

■¥kin = a i { X )x l , (4.14) 
where we introduced, following the notations in [7], the auxiliary notation for the fermonic fields: 



XI = i>l , X2 = i>2 , X3 = Ipl , Xi = i>l ■ 
In this case, the functions dj(x) take the form: 



a i 



a 2 



«3 



04 



i iVXg 2 , iVXg3 , 

2 XS ZP~ X3XiXl + gj3 X1X2X3X4.X4, 

i i\f\g2 , . iVXg3 , 

2 X4 ^2^^3X4X2 + 8J3 XlX2X3X4X3> 

i\f\g 2 , iV\gs , 

Xi 4 2J2 X1X2X3 + 8J3 X1X2X3X4X21 



1 
2 
i 

2 X2+ 2 J 2 



92 



X1X2X4 + 



93 



8J 3 



-X1X2X3X4X1- 



(4.15) 

(4.16) 
(4.17) 
(4.18) 
(4.19) 



The next step is to derive the Euler-Lagrange equations following from (4.12). Thus, we consider: 

5(J?) = 5[a l ( X )xi]-H^) = 0, (4.20) 

with the implied sum over i going from 1 to 4. 

It is easier to evaluate each variation separately: 



S[ai (x)xi] = / dy{ s Xj(y) 

5aj(x 



Sajjx) _ 6aj(x) 

$Xj(y) v Sx'j{y) 



Xi(x) 



+ $Xj(x) 



< , sXi{y) + pjr\ Xi{y) 
oxi{y) SXjiy) 



(4.21) 



To make it possible to write the Euler-Lagrange equations as in (4.8), we must be able to write this 
variation as follows: 

5 [oj (x) Xi] = $Xj(x) fy^x) Xi{x). (4.22) 

Clearly, the expression (4.21) does not have such a form. However, as we will show bellow, once we 
fix the values of the indices i and j, it is possible to reduce (4.21) to (4.22). Let us then compute 
the right hand side of (4.21) for the case i = j = 1. One has: 



dy { <5xi(y) 



Sxi(x) 



Sai(x) _ Sai(x) 

$xi{y) y $x'i(y) 
iV\g2 



xi(x) + 5xi 0*0 



Sa^x) . 6at(x) ., 

6^) Xl{y) + 5xW) Xl{y \ 



2J 2 

5xi(x) n n (x) xi(x) 



WXg3 
4J 3 



(X3X4 - X4XQ + A T 3 X2X3X4X4 



Xl{X) 



with 



^11 0*0 = - ^~ 1X3X4 - X4X 3 j + A T 3 X2X3X4X4- 



, , : i VX g 3 
4J 3 



(4.23) 
(4.24) 



WX ff2 
2J 2 

Repeating this calculation for the other indices we can obtain all the elements of the matrix 

Qij(x), which we collect in appendix E. 



15 



Computing the variation of Jrif one obtains: 



5jr(x) = / dy 



Sxi(y) -t^t + ^ y ir>H 
oxi{y) Wily) 



S Xi {x) Hi{x), (4.25) 



where we have introduced the functions: 



Hi = —X4~X3, H 2 = — J-X3 + X4, H 3 = — J-X2 + X1 and H 4 = — Xi ~ Xi- 

Now, substituting (4.22) and (4.25) back into the equation (4.20), we obtain: 

n ij (x) Xj (x) = H i (x), (4.26) 

which is in direct correspondence to (4.8). 

The Dirac structure is defined in the standard manner [3]. Let F [ X i( x )] an d G [Xi{ x )\ be two 
functionals of the fields and define the Dirac brackets between them in the usual way: 

K fl ) = l""^»»WW (4 ' 27) 

where oj 13 is some function of x, y G R. For any even functional of the fields F [xi( x )]i we can write: 

r 5F 

{ X k(z),F} = / dw { Xk (z), X i(w)} (4.28) 
J oxi{w) 

Clearly, the Hamiltonian (4.5) is an even functional of the fields, so that we can use (4.28) to write 
Hamilton's equations as: 



Xi( x ) = {Xi(x),H} = J dy {Xi( x ),^(y)} = Jj dy dz {xi( x ),Xj(z)} Hj(y)5(y - z) 

dy{xi{x) lX 3{y)}Hj{y)- (4-29) 



Finally, since the matrix Q(x) is non-singular, 11 one can invert it and combine the equations (4.26) 
and (4.29) as follows: 



Xj {x) = nrl(x)Hi(x) = J dy {Xi( x ),Xj(y)} ^{y). (4.30) 

Then, we obtain: 

{Xi(x),Xj(y)} = tl£(x)6(x - y). (4.31) 

Thus, the matrix f2 _1 (a;) defines the Dirac structure. We collect all the matrix elements of Vt^ 1 {x) 
in appendix F. 



11 The non-singularity of Q(x) can be directly established by explicitly checking the existence of the inverse matrix 
fl^ 1 (x). See appendices E and F for details. 
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5 The Algebra of Lax operators 



With the reduced 2x2 Lax connection derived in section 3 and the Dirac brackets deduced in 
section 4, we are finally in the position to obtain the algebra between the Lax operators and show 
its non-ultralocal structure. We start by also rescaling the fields in the expression for the spacial 
component of the Lax connection (3.1): \i ~^ J^ 2 Xii so that it is consistent with the Lagrangian 
(4.1), which we used to derive the Dirac algebra, provided we further impose the constraint (2.2) 
by setting 52 = 1 = 93- 12 Moreover, we decompose the spacial component of the Lax connection in 
a more convenient structure: 



£i(x;m) = ^ a \x; fi)I + ^[ a \x; fi)J + A a (x; 

= 4 CT) (x; /i)l 2 + ^ (x; fi)a 3 + AH ( x; ^ + A« (x; fi) 



(7 



(5.1) 



where a 1 , i = +, — , 3 correspond to the usual Pauli matrices. Here we also introduced the functions 
£j a \x; (J,), j = 0, 1 and A^\x; ll). The former (latter) are even (odd) polynomials of the fermionic 
fields, containing at most one space-derivative, the expressions of which are relegated to appendix 
G. 

Using the decomposition (5.1), we reduce the task of computing the Dirac brackets between two 
L-operators to the evaluation of the following sixteen Dirac brackets between the functions £j°^(x; ll) 
and Act (a;; ll): 



{£i(*;mi) ? = {&\x-» 1 )^\y-Li 2 )) 1 2 ® 1 2 + [&\x; Lix)^\y; ^2)} 1 2 

{^ CT) (x; W ), A(-)(y;/i 2 )} 1 2 ® a + + {^ ct) (x; Ml), A«(y; /x 2 )} 1 2 ® o~ 

{cf'Wi). ^ M2)} a 3 ® i 2 + {ei ff) (x; mi), ^(y; /x 2 )} a 3 ® a 3 

{^ ) (x;/, 1 ),A(-)(y; / , 2 )} ( T 3 » C T + + {e i ,7) (x; / u 1 ),AW(y; / , 2 )}a 3 »a- 

{ AH (x; m), (y; M2) } <x + ® ^ + + { ( x; A^ ( y; M2 ) J a + g, 
{AW Or; Mi), ^ fa M2) } ^~ ® 1 2 + {AW (x; M i), ^ (v, M2) } <r~ ® a 3 
{ A (+) (x; mi), AH (y; M2) } <r _ ^ + + { A « (x; Ml), A^ (y; M2) } a" ® «7". 



(7 



+ 
+ 
+ 
+ 
+ 
+ 
+ 



(5.2) 



At first glance it seems that this decomposition only makes the calculations even more daunting. 
However, this computation is, in fact, easier, since only half of the Dirac brackets need to be eval- 
uated. The other brackets can be obtained by taking the involution of the corresponding brackets, 
as we will explain bellow. 

Let A (xi(x)) and B (xi(y)) be two arbitrary functions of the fields, then the behaviour of the 
Dirac brackets defined by the matrix O -1 through (4.31) under involution is determined by the 
parity of such functions, 



{A(x),B(y)Y 



+ {A*(x),B*(y)} , if A or B even , 
- {A* (x) , B* (y) } , if A and B odd . 



(5.3) 



12 Alternatively, we could use the inverse scaling transformations (see footnote 3) to restore the general coupling 
constants gi and gs satisfying the constraint (2.2). 
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Furthermore, taking into account that the functions Q (x; /i) and (x; /i) have a defined parity 
and behave under involution as: 

i i f ) \x^) = -i ( f\x^*) and AH*(x, M )=AW(x, M *), (5.4) 

one can obtain some non-trivial relations amongst the following brackets: 

{^\x^ 1 ),A^\y, f , 2 )Y = -{^\x-,^),A^(y,^ 2 )}, (5.5) 
{4 ± )(x,/ il ),AW(y,^)}* = {AW(x,^),A^)(y,^)}. (5.6) 



Moreover, if we restrict the arbitrary functions of the fields ^4 (xi(x)) and B (Xi(y)) to the subset 

:(<■■ 
j 



comprised of the functions Q (x; /i) and A^ (x; /x), one can show that the Dirac brackets between 



A (xi(x)) and I? (Xi{y)) have the general form: 13 

{A(x),B(y)} = h(x)5(x - y) + f 2 (x)d s 6(x - y) + h(x)d 2 J{x - y), (5.7) 

where fi(x), i = 1,2,3 are some polynomials of the fields and their space derivatives. Then, by 
invoking the (anti)symmetry of the brackets: 

1 V h yyn \ + {B(y),A{x)} , if A and B odd , V ' 

one can easily derive the following relation: 

{B(x), A(y)} = T { [fl(x) - d x f 2 (x) + d 2 x h(x)] S(x - y) + [-/ 2 (x) + 2d x f 3 (x)] d x S(x - y) 

+ f 3 ( x )d 2 J(x - y)} , (5.9) 

where the /j(x) are the same functions appearing in (5.7). Here the plus sign corresponds to the 
case where both fields are odd and the minus sign to the case where at least one of the fields is even. 
Hence, even though there is no simple equation that relates {A(x),B(y)} to {B(x), A(y)}, we have 
showed that by knowing the former one can readily obtain the expression for the latter. 

By using the properties described by the equations (5.5) and (5.6), and by the relations (5.7) 
and (5.9), we can drastically reduce the amount of calculations necessary to derive the algebra of 
the Lax operators. Namely, out of the sixteen brackets in equation (5.2), only seven are actually 
independent. Besides that, these properties are also extremely useful for computing each of these 
independent brackets. 

Before proceeding with the calculation, it is worth noting that both (5.7) and (5.9) already 
display a severe non-ultralocality, as both of them, not only contain the first derivative of the delta 
function, but also its second derivative. Given the structure of (5.2), one could still hope that 
some non-trivial cancelations would still render the algebra ultralocal, or at least diminish its non- 
ultralocality. However, this is not the case for the AAF model, and, as we will show bellow, they 
remain even after all the terms are added together. 

13 The relevant question one could still pose about the generality of the expression (5.7) regards the truncation of 
this series at the second derivative of the delta function. As a matter of fact, there cannot be terms proportional to 
the third derivative of the delta function or higher, because the functions ^ CT ' [x; /i) and {x; /i) are at most linear 
in the space derivatives. Hence, at most two space derivatives can act on the delta function after the Dirac brackets 
are computed. We note, nevertheless, that the functions fi(x) can still be proportional to some derivative of the 
fields. In this case the derivatives may come not only from the functions Q (x; /«) and A^ (x; fj), but also from the 
Dirac structure itself. 
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Computing all the independent brackets in (5.2) is still a very lengthy and tedious task. We 
sketch here only the main steps of this computation for one of the independent brackets. For the 
sake of clarity, let us consider the first brackets in (5.2). We first decompose it in a sum of Dirac 
brackets such that each entry of the Dirac brackets contains only monomials of the fields at some 
given order. For instance, 



j2 {A\(v\V-i),Ax(y;n2)} + {Ai(x; m), A 2 {y; ^ 2 )} 
+ jg {A 2 {x; m), A^y; fi 2 )} + {A 2 (x; /j, 1 ),A 2 {y; fi 2 )} 



(5.10) 



where A\ contains only quadratic terms in the fields and its derivatives, and A 2 only quartic: 

A\ = -X3Xi + X4X2 ~ X1X3 + X2X4, (5-11) 
A 2 = X2X3X4X1 + X1X3X4X2 + X1X2X4X3 - X1X2X3X4- ( 5 - 12 ) 



It is easy, though tiresome, to compute a complete set of all the Dirac brackets between all 
possible combinations of monomials of the fields and its derivatives. For example, the one needed 
to compute {Ai(x; Hi), A\{y\ ^2)} is: 



{xi{x)xj{x),xk{y)x[{y)} 



d x 



XiX'k^jl ~ XiXfojk 



- XiXk^ji 1 - XiXi^jk 1 - XjXk^ii 1 - XjXi^ik} S(x - y) 

+ {d x [xiXk^]f\ + [Xkx'i + XiXfc] fij, 1 + XiX'i^t + XkX'Pu} 9 * 6 ( x ~ y) 



+ XiXkM jl 1 d x 5(x - y). 
After multiple substitutions and summations one obtains the expressions of the type: 



{A 1 {x;n 1 ),A x (y;ii 2 )} 
{Ax(x;ni),A 2 {y;n2)} 



16 
l r 



[2d x T n 6(x -y)+ WndJix - y)] , 



T$5(x-y) + T[1>d x 5(x-y) 



(2). 



(5.13) 

(5.14) 
(5.15) 



where the explicit expressions for Tn, and are presented in appendix G. 

Deriving the expression for {A 2 (x; fj,x), A\(y; fi 2 )} becomes a straightforward task by employing 
the property displayed by equations (5.7) and (5.9), since we can read off the concrete formulas for 
the fi(x) by comparing (5.7) to (5.15) and then simply substitute them into (5.9) to obtain: 

1 



{A 2 {x\m),Ax(y;n 2 )} = - { 



(2) 
x± 12 



6(x-y) + r ( $d x 6(x-y)}. (5.16) 



Finally, we consider the last brackets in (5.10). This calculation becomes rather simple, by realizing 
that the relevant Dirac brackets have the form: 



{xi(x)xj(x)xi{x)xrn(x),Xn(y)xp(y)xq(y)x'r(.y)} = h{x)S(x-y). 



(5.17) 



Thus, 



{A 2 (x; m),A 2 (y; /x 2 )} = F(x)5(x - y), (5.18) 

where both fi(x) and F(x) are some polynomials in the fields and their spacial derivatives. Since 
A 2 is an even function of the fermionic fields, we can use property outlined by the equations (5.7) 
and (5.9) to write: 

{A 2 (y; »i),A 2 (x; m )} = -F(x)5(x - y). (5.19) 
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On the other hand, we can use the fact that A 2 does not depend on the spectral parameter, and 
one can simply exchange x -H- y in the equation (5.18) to conclude that: 



{A 2 {x;ti 1 ),A 2 (y;fi2)} = 0. 



(5.20) 



It is interesting to stress that the explicit form of the matrix elements O^- 1 was not used in the 



derivation of (5.14), (5.15), (5.16) nor (5.20), that is, Tn, T^ anci ^12 are ST,m wr itt en in terms of 



,(2) 



This means that the non-ultralocal form of the decompositions (5.14), (5.15) and (5.16) are 

not a direct consequence of the Dirac structure of the model, but, up to this point, only of the form 

(2) 

of the Lax connection. Thus, were Tu and T\ 2 to vanish identically upon the substitutions of the 
expressions for fiL 1 , the equations (5.14), (5.15) and (5.16) would be ultralocal. However, this is 
not the case. In fact, after careful substitution of the expressions for f^. 1 , given in appendix F, we 
conclude that: 



16J 2 



2T n + 



.(2) 
12 

J 



6(x -y) + 2 



2T 



11 



^(2) 
■ 12 

J 



d x 5(x - y) 



(5.21) 
1 (2) 

where the sum 2Fn + jT\ 2 extends up to the sixth order in the fermion and its spacial derivatives. 
We relegate this lengthy expression to appendix G. 

The remaining Dirac brackets in (5.2) can be computed in exactly the same vein. In fact, the 
brackets which we have just computed are one of the simplest. Some of them are even plagued by 
more severe non-ultralocality, which extends up to the second derivative of the delta function, in 
full consonance with what one would expect from (5.7). It turns out that even after summing all 
the brackets according to (5.2), this severe non-ultralocality survives. In particular, this prevents 
one from using the standard methods in the context of integrable models, as it is not possible, in 
general, to introduce a well-defined r-matrix formulation starting from a non-ultralocal algebra. 11 
Nevertheless, there exists a general framework to deal with such non-ultralocal systems, which was 
introduced by Maillet in [22-25], 15 generalization of which, suitable for the ^4^4.F-model, we consider 
in section 6. 



5.1 Fermionic Wadati Model 

Although the AAF model provides a very representative and interesting example of highly non- 
ultralocal models, its algebra is by far too complicated for an initial analysis. Therefore, it is more 
illuminating to consider first a simpler model. One such model can be obtained by considering the 
field independent truncation of the AAF Dirac algebra, which, as we show below, corresponds to 
the fermionic version of the Wadati model [41,42]. In this case, the intricate Dirac structure O^ 1 
presented in appendix F simplifies considerably, and the resulting Dirac brackets are the canonical 
ones: 

{Xi(x),X3(y)} =iS(x-y), , 522 , 
{X2(x),X4(y)} = iS(x -y). 

However, despite the simplicity of the Dirac brackets, the algebra between the Lax operators 
(5.2) is still quite complicated, as a result of the form of the Lax connection itself. 16 Substitution 

14 There are, however, some exceptions for which the non-ultralocal algebra for //-operators still leads to a well- 
defined algebra for the monodromy matrices. See the footnote 2. 
15 See also the earlier attempts by Tsyplaev [41]. 
16 See the discussion bellow (5.20). 
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of the canonical Dirac structure (5.22) into the independent brackets in (5.2) corresponds to their 
truncation to the first order in the fermion: 



\i a) (x; m) , Ai +) (y; n 2 )\ = —3 [-^hi^xi + k^Xs] d x 5(x - y) 

[ih{^2)XA ~ U{P2)X?>] dld{x - y), 
1 / 2J 



(5.23) 



1 

+ — 

4J2 



{^\x^ 1 ),Ai + Hy;^)} 



8Jf \VX 

[h(fj,i)h(fJ>2) + hi^hi^)] X3 ~ i [^1(^1)^3(^2) + ^2(^1)^4(^2)] Xa 9 x S(x - y) 



(5.24) 



H — ~ 3* ( [h{ni)h(fJ-2) + ^2(^1)^3(^2)] X3 + i [^1(^1)^3(^2) + ^2(^1)^4(^2)] Xa ) dlK x ~ 2/)> 



{AS >(x;//i), A£%;/x 2 )} = -^M^i)h{^2) + M^i)U{^2)]dl5{x - y). 



(5.25) 



The remaining brackets from (5.2) are either identically zero at this order, or can be trivially derived 
through the use of properties (5.5), (5.6), (5.7) and (5.9). 

Collecting everything together, we can write the algebra between the Lax operators as follows: 



{Ci(x;m) f A(y;/x 2 )} = N (x, y; m, fM 2 )5(x - y) + Ni(x,y; /J, 2 )d x 5(x - y) 

+ N 2 {x,y;fi 1 ,fi 2 )dl5(x - y), 



with the matrices: 



N (x,y;iii,ii 2 ) 



Ni(x,y;m,ii 2 ) 



N 2 (x,y;fj ll ,fj, 2 ) 






A^Vi,/^) 




^Vi (1) (w,M2) 
^1 (2) (W,M2) 




AfVl,M2) 
r(2) 









-Jv (1) *(mi,/4) 





(2)* 





-iV (2) *Gu|,|U* 2 ) 














M 3) (ti,i4) 




-iV 2 (3) *(^,/i|) 
-iV 2 (4) *(^, M |) 





AfVi,^) ^! (4) (/xi,M2) 



AT«*(^,/4) -iV^ (M*i,/4) 

^2 (5) (M1,M2) 



( 4 )/ 



(2)*. 







V 







JV 2 (3) ( MllM2 ) iY^( m ,/i 2 ) 



The functions (fj,i, fi 2 ) are defined in appendix H. 

Thus, restricting to the field-independent part in the algebra (5.26), we find: 

{£i(x; A) ? City] P2)} = -j cosh (ft + /3 2 ) [<r+ ® a~ + a" ® <r+] <9^(x - y), 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



(5.30) 



where we have used the parametrization (A. 5) for the U functions. The algebra (5.30) has, curiously, 
the same structure as the Wadati model, and therefore, can be considered its fermionic counterpart. 
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We refer to the original literature [41,42] for more details and discussions of the Wadati model. 
We only mention here, that despite the highly non-ultralocal algebra, the classical r-matrix in the 
infinite space limit, nevertheless, can be found. It is given, for example, in [41]. Although the full 
algebra for the AAF-model is considerably more complicated than the truncated algebra (5.30), it 
is still an interesting example to consider, as it demonstrates the necessity of new methods without 
technical complications. 

6 Generalized Maillet algebra 

In this section we will briefly explain the Maillet 's formalism of the r- and s-matrices, and the 
symmetric limit prescription to deal with the non-ultralocalities in the algebra. Afterwards we shall 
develop its generalization suitable for the ^^.F-model. The latter is needed, since, as we have seen 
in the previous section, the algebra for the AAF model contains terms up to the second derivative 
of the delta function. The more complete account of such generalized Maillet algebras will be given 
in the separate publication [43]. 

The starting point of the formalism considered by Maillet is the non-ultralocal algebra of the 
form: 1 

{C(z; \)®£(z'; ft)} = A{z, z'; \)8{z - z') + B(z, z'; X)d z >5(z - z') + C(z, z'; X)d z 5(z - z'), (6.1) 

which can be recast in the following more convenient form: 

{C{z; A)®£(/; ft)} = (d z r(z; A, ft) + [r(z; A, ft),£(z; A) 1 + 10 C(z; ft)] 
+ [s(z; A, ft), 1 £{z; ft) - C(z; A) 1]) 5(z - z ) 
+ s(z;X,ft)d zl 6(z- z') + s(z';X,fi)d z 6(z- z'). (6.2) 

The r- and the new s-matrices are defined as follows: 

s(z; A, ft) = - [B(z, z; A, ft) - C(z, z; A, fi)] , (6.3) 
r(z; A, fi) = - [B(z, z; A, ft) + C{z, z; A, ft)] + r (z; A, ft), (6.4) 

with the function tq(z; X, ft) in (6.4) being determined from the equation: 

d z r (z; A, ft) + [r (z; A, ft),£(z; A) <g> 1 + 10 C(z, ft)] = Q(z; A, ft), (6.5) 

where the function Q(z; A, ft) has the form: 

Q(z; A, ft) = A(z; A, ft) - d v [B(z, v; A, ft) + C(v, z; X, ft)] v=z 

+ [10 C(x; ft),B(z, z; A, ft)] + [£(x; A) 1, C(z, z; X, ft)] . (6.6) 

Thus, the algebra (6.2) is defined now by the two matrices r(z; A, ft) and s(z; A, ft), which can now, 
in general, depend on the coordinates. 

17 The splitting into B(z, z 1 '; A, n) and C(z, «'; A, fx) functions is introduced for the sake of convenience. One may 
choose these functions (as it was done for example in the original paper [22]), so that PB(z, z';X,fi)P = — C(z' , z; /x, A), 
where P is the permutation operator. This, however, is not necessary, and the choice can be a different one. 
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Using these definitions, it is easy to derive the algebra for the transition matrices T(x,x'; A): 



{T(x,y;\)®T(x',y';»)} 

= T(x, x ; A) (g> T(x\ x ; li) r(x ; a > A 4 ) T (xo, y, A) &> T(x , y; £t) 

- ^(a;, y ; A) 8) T(x', y ; n) r(y ; A, li) T(y , 2/; A) <g> T(y ,y'; li) 

+ e(x - x')T(x, x ; A) (g) T(x', x ; /u) s(x ; A, /i) T(x , y; A) (8) T(x , y; n) 

- e(y - v)T{x, y ; A) ® T(x, y ; n) s(y ; A, /u) T(y ,y; A) (g> T(y , 2/'; li), (6.7) 

where all the x,x',y,y' are distinct, xo = min(x,x'); y$ = max(y,y'), and the ordering is chosen 
such that x, x' > y,y'- The crucial point in this formula is the appearance of the functions e(x — x') 
and e(y — y'), which make the algebra (6.7) ill-defined. Indeed, the discontinuity in the algebra (6.7) 
at coinciding points x = x' (or y = y') is proportional to the value of the s-matrix at this point. 
Furthermore, as shown in [22], it is not possible to define the Poisson brackets at the coinciding 
points, such that the Jacobi identity is satisfied. However, it is possible to define weak Poisson 
brackets, via the Maillet's symmetric limit procedure, such that the Jacobi identity, as well as all 
the standard properties hold. 

The symmetrization procedure and the weak Poisson brackets are defined for each n-nested 
brackets of the type: 

A n { Xi , y t - \) := {T(x lt yi\ Ai)f {• • • f {T(x n , y n ; A n )® T(x n+1 , y n+1 ; A n )} ...}}. (6.8) 

This quantity is only well-defined if all the points X{ and yi are distinct. For coinciding points, 
one introduces the weak Maillet brackets by a point-splitting and symmetrization procedure. For 
example, for X{ = x, one defines: 

A n (x, yi ; Xi) := lim - V A n (x + ecr(l), . . . , x + ea(n + l), yi ; A;) , (6.9) 
e->0 in + 1) 

where P stands for all possible permutations of (1, . . . , n + 1). In particular, for n = 2, this definition 
leads to the Maillet brackets between the transition matrices: 

{T(£,y;A)®T(x,7/ ;Ai )}M 
:= \ lim ({TO - e, y; A)? T(x + e, y'; fx)} + {T{x + e, y; A)f T{x - e, y'; fj,)}) . (6.10) 

With such well-defined weak brackets, one can show that the algebras for the transition matrices 
for the equal and adjacent intervals have the form: 

{T(x, y; A)f T(x, y; n)} M = r(x; A, fj,) T(x, y; A) (g) T(x, y; fj,) (6.11) 

- T(x, y; A) 8) T(x, y; fx) r(y; A, li), 
{T(x, y; A)f T(y, z; li)} m = {T(x, y; A) ® 1) s{y; A, li) (1 ® T(y, z; A)). (6.12) 

We refer the reader to the original papers [22-25] for detailed discussion and several impressive 
examples, which can be solved with this prescription. Moreover, it is possible to formulate these 
relations on the lattice, which makes it possible, in principle, to proceed with quantization [44]. 

6.1 Operator regularization method and the classical limit 

Before considering the generalization of these ideas for the AAF model, we address the validity 
of using such weakly defined Poisson brackets in the classical theory. Indeed, the construction of 
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the Maillet brackets via the symmetric limit prescription is not intuitively plausible, and seems 
somewhat artificial. We claim, however, that Maillet's prescription is the natural consequence of 
the regularized quantum theory, using the split-point regularization introduced in [19,20]. Although 
it was demonstrated on the particular example of the Landau-Lifshitz model, the method itself is 
general enough to be applicable to other singular systems, and, therefore, we consider here only the 
general ideas, and leave the complete proof for the future publication [43]. 

The central concept in the operator regularization method is the introduction of the regularized 
operators Sjr(x): 

S a Ax):= J dZT e (x,0S a (0, (6.13) 

where J r e ^x,^) is some symmetric and smooth function, which depends on some parameter e. The 
function J- e (x, £) should also satisfy some additional conditions depending on the algebraic structure 
of the model, such that the singular expressions (for example, the Yang-Baxter equation) due to 
the operators product at the same point, become well-defined. Moreover, in the e — > one should 
restore the original expressions. The theory can then be reformulated in terms of the regularized 
operators Sjr(x), and only in the end one should remove the regularization e — > 0. This program 
has been completed for the Landau-Lifshitz model, and it was shown that the spectrum and the 
higher order charges can be obtained by using this method [19,20]. In addition, all the quantum 
charges were shown to be self-adjoint, and the necessary self-adjoint extensions were constructed. 
We also established a connection between the self-adointness of the quantum charges and the S- 
matrix factorization. It is clear, that the function J- e (x,tf) in (6.13) is essentially equivalent to the 
split-point procedure in the quantum theory, so that the algebra, and other constraints become well- 
defined. One should consider the quantum theory as fundamental, and the classical theory should 
follow from the quasi-classical limit: H — )■ 0. Therefore, such split-point operator regularization in 
the quantum theory, which makes all the singular expressions well defined, should also appear in 
the quasi-classical limit. The Maillet symmetric limit prescription (6.9) is exactly such split-point 
regularization, and we conclude that the Maillet's definition of the Poisson brackets is simply the 
consequence of the regularized operators of the quantum theory. The full proof and details of this 
argument will be presented in a future publication. 

6.2 Maillet algebra for the AAF model 

We turn now to the algebra for the L-operators for the AAF model. As we have seen in the 
previous section, the algebra has a more general form, and contains terms up to the second derivative 
of the delta function. Namely, the algebra has the following form: 

2 

{C(z- \)®C{z'; fi)} = A(z, z'; X)5(z - z') + B VJ {z, z'; A)d*dJ,<5(z - z'). (6.14) 

ij=0 

Although we have found the A and £?« functions, we will not present here their explicit expressions, 
due to their lengthy form. Instead, we will focus on the some of the preliminary essential conse- 
quences of the algebraic structure (6.14) for the bosonic case with arbitrary functions A and Bij, 
leaving the full details and explicit expressions to a separate publication. 
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One can recast the algebra (6.14) in a more convenient form, similar to the one in (6.2): 



{C{z- \)®C{z'- fj)} = (d z r(z; A, y) + [r(z; A, //), C(z; A) ® 1 + 10 £(z; /x)] 
+ [ai(z; A, /i), 1 (8) C(z\ n) - C(z; A) (g> 1] 
+ [d z s 2 (z; A, ft), £(«; A) (8) 1 + 10 £(z; //)] 
+ [[* 2 (*; A, a*), A) 1] , 1 /x)] 
+ [[aaCs; A, //), 1 ® £(z; /j)] , C(z; A) ® 1]) S(z - z) 
+ si(z; A, fj,)d z /5(z — z') + s\{z'; A, y)d z 5(z — z') 
+ s 2 (z; \,y)d 2 z ,5(z - z) + s 2 (z'; \,y)d 2 z 5(z - z), (6.15) 

where the r, si and s 2 matrices are defined using the functions A,Bij, similar to (6.3)-(6.6). The 
essential difference is the introduction of the additional matrix s 2 , due to the second derivative term 
in the algebra (6.14). It is easy to show that the algebra for the transition matrices, that follows 
from (6.15) has the form: 

{T(x,y;X)®T(x',y';iJ,)} 

= T(x, x ; A) (g> T(x', x ; fj,) u(x ; A, fx) T(x , y; A) T(x , y'\ y) 

- T(x, y ; A) T(x', y ; //) u(y ; A, ft) T(y , y; A) (g> T(y ,y'; /i) 

+ e(x - x')T(x, x ; A) T(x', x ; /i) u (x ; A, /j) T(x , y; A) ® T(x , y'; /j) 

- e(y' - y)T(s, y ; A) T(x', y ; m) ^(yo; A, /i) T(y , y; A) r(y , y'; /j), (6.16) 

where we have defined: 

u(x; A, /i) = r(x; A, /i) + 9 s s 2 (x; A, /j) + [s 2 (x; A, /i), A) ® 1 + 1 £(x; /j)] (6.17) 
u(x; A, /i) = si(x; A, /i) + [s 2 (x; A, [i),£(x; A) <g> 1 — 1 ® £(x; //)] (6.18) 

This algebra has exactly the same form as the algebra (6.7), where the r(x; A,/i) and s(x; A, ft) 
functions are replaced by the u(x; A, /i) and v(x; A, /i) functions. Hence, one can immediately write 
the algebra for transition matrices for the equal and adjacent intervals from (6.16) by replacing r — > u 
and s — > v in (6.11) and (6.12), as well as use the general results of the previous sections, including 
the construction of the symmetric limit prescription for nested brackets. Thus, we conclude that 
even though the algebra for L-operators (6.15) is modified by the additional terms depending on 
the s 2 (x; A,/i) matrix, the resulting algebra for the transition matrices (6.16) is the same, and the 
complication introduced by the second derivative term is only technical. The essential steps of the 
Maillet prescription can be repeated without any conceptual complications for the graded algebra 
as well, which is the case of the AAF model, the generalization is straightforward. Finally, we note 
that the lattice version of the algebra (6.16) is known as well [44], and may in principle be used to 
construct the lattice version of the AAF model. The full details of this analysis will be presented 
in the future publication [43]. 

7 Conclusion 

We have considered in this paper the classical integrability of the AAF model, starting from the 
first principles of the inverse scattering method. The necessity of this consideration was stimulated 
by the fact that the perturbative calculation still do not provide the complete information about the 
integrable properties of the AAF model. For example, the perturbative approach was not sensitive 
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to the rather complicated structure of the Dirac brackets, which, therefore, was not taken into 
account. In the process of our consideration, we have proved that the Lax operator surprisingly 
admits a simpler 2x2 representation. This is especially curious, since the much simpler Thirring 
model has a more complex 3x3 known representation. We have also found that the constraint on 
the coupling constants g\ = <?3, which was derived in [14] from the S'-matrix factorization condition, 
must also be satisfied in the classical theory, and that for no extension of the Lax connections 
integr ability holds without this constraint. 

The algebraic structure of the AAF model, however, is rather complicated and has a highly 
non-ultralocal form, which contains terms up to the second derivative in delta function in the 
algebra of the Lax operators. This led us to the generalization of the regularization technique due 
to Maillet. Here, it was necessary to introduce three matrices r, s\ and S2 to properly encode the 
algebraic structure. We have also derived the algebra for transition matrices for the AAF model, 
and have shown it to have exactly the same form as it is the case of the usual non-ultralocality, with 
modified r- and s-matrices. Strictly speaking, although we have obtained all the general expressions, 
further progress depends on the analysis of the equation defining the r-matrix, and finding its local 
solutions. This will allow one to proceed in the manner similar to, for example, the complex sine- 
Gordon model, where the r-matrix has indeed a local, coordinate independent form, and the model 
can be solved exactly in the infinite space limit [22]. 

We also obtained in the process the fermionic counterpart of the Wadati model, which introduces 
an interesting toy-model for testing the generalized Maillet algebra, as it seems possible to explicitly 
compute some relevant quantities, such as the transition matrix and its infinite line limit. 

We have also proposed that the symmetric limit prescription of Maillet has its origin in the 
regularized quantum theory. Indeed, based on the method proposed in [20], the quantum theory is 
made well-defined via operator regularization. Maillet algebra is then obtained in the quasi-classical 
limit. This is an interesting direction which will be explored in more details in the future. 

The problem of obtaining an alternative Lax connection, for which the algebra becomes ultralocal 
is open. In principle there may exist such a connection, which does not follow from the reduction 
of the strings on AdS§ x S 5 to the su(l, 1) subsector. It would also be interesting to obtain the 
equivalent model via bosonization procedure. Indeed, it is easy to see, by examining the constraints 
of the gauge fixing procedure in [7], that the additional higher order non-ultralocality appears in 
the algebra due to eliminating the bosonic fields in favor of the fermionic ones. Thus, it is very 
desirable to obtain the bosonic version of the AAF model, which should be less non-ultralocal, and 
may make the integrability analysis a much simpler task. 

Appendices 
A Notations 

We use the following representation for the two-dimensional Dirac matrices: 





Following the notations in [7], we choose here the sign of k as: 




(A.2) 



K = 



2 
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In this case the functions l{ have the following dependence on the spectral parameter \x (see the 
discussion in [7, 12]): 



Iq — 1, l\ 



, h = s\ 



2// 



1 — fl z \ — ^ 

with the following choice of parameters: 



, h = $2 



and 1 4 = S3 



(A.3) 



(A.4) 



si + S2S3 = and (S2) 2 = (S3) 2 = 1- 
Alternatively, one could use the parametrization: 

Z = l, h = cosh(2/3), Z 2 = -2tanh(/3), i 3 = cosh(/9) and Z 4 = sinh(/3). (A. 5) 

B Equations of Motion in the reduced form 

In this appendix we give the equations of motion for the AAF model in the reduced form. As 
explained in the main text, they are derived from the original equations of motion (2.3) and (2.3) 
by recursively eliminating the time derivatives present in the cubic and higher order terms, through 
multiple usage of the equations of motion. 

The reduced equation for ipi, following from the AAF Lagrangian, has the form: 



iJipi — VXtp'2 + Jipi 



A52 



2J 



+J (tpiip&i - ^2^2 i>i + iPi^'i - V'iV'iV^)] 



K92? 

4J 2 



- ^ {^l^l^l^'l ~ 3^2*^2 ^1^2 + 2 V , lV'2V ; 1^2V ; I - ^\i>\^\^2^2 

+^2*V>l^l) +2J^ 2 VlV'2^2] + (A!> ~ {mY ' 

(X) 3/2 92g 3 



4J 2 



8J 3 



8V A^> 2 ^ ifeV'l "02^2^1 

0. 



(B.l) 



while the reduced equation for ip2i which follows from the AAF Lagrangian, is: 

%2 



1J1P2 - VXtp'x + Jlp2 



2,1 



Vx (-2^1 + ^2^2 + - ^TM'i + 

A (S^iriM 



K92? 



+ J {tfefa^l - ^2 ^2^2 - ^2^1 ^2 + ^l^iVl)] + 4J2 
+3V4Vl^2^i + 2^'*^2^t + ^2*^2*^2 + V4< Wl^) + J (-3^2*^2^ 



8J 



A (3^f v^iV^'i + 3*0^2 ^\i>2$2 - 2^1 V^V^ + V4* ^\i>\^%i>\ 

(A) 3 / 2 ( 52 ) 3 



8J 3 



4J 2 



0. 



(B.2) 
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The dynamical equation for 0i, following from the zero-curvature condition, has the form: 



iJipi — V~\tp2 + Jipi 

-020101 + 02*0101 
'A 



4 

U 



+ 



2020202 - 010102 + 201*0102 ~ 0*0102 ~ 010102 + 01*0102 

A 



20|0 2 01 + 010101 - 02 0201 

3J r , , , A 



4J 



[020201 - 01*0101 ~ 02* 0201 ] 



-4V>202010102 - 0101020202 



+ "J" [020101 + 020202 + 01* 0201 J ~ ^ [02020.., | (j y 
+0202010102 + 0101020202 - 01*01'020202 ~ 02*020*0102 + 2lp*lj/ 1 1p*1pli>2 ~ 0* 01 01*0102 
+2^*01020202 + 2 *'010 2 0201 + 0201010101 ' 



16J 



[-30*;020*010i - 01*01020201 



+02*0201*0101 " 20 2 *0 2 0(*0 2 0i] = 0, 

while the dynamical equation for 2 , which follows from the zero-curvature condition, is: 

iVA 



(B.3) 



i J02 - n/Ai/4 + </02 



-0*0 2 02 + 01*0202 



+ 



iJ 



2ll){l}) 1 1p' 1 - 02020'l + 20 2 *0 2 01 - 020201 - 020201 + 02*0201 

A 



2^1 01 ^2 + 020202 - ^10102 



4J 



[010102 - 02*0202 - 01*0102] 



3J 



[010202 + 02*0102 + 01010i] + -p [01010: 



A 



40i0i02020i + 2^202010101 



-010102020! + 0202010101 + 0101020201 + 02*01010102 + 0i* 0102 0201 ~ 2020 2 020201 
+02*02020201 - 0102020202 ~ 0i*02010102 1 



16J 



[-30i0'i020202 - 02*02010102 



+01*0X0^0202 " 2V>i0i01*0102] 



0. 



(B.4) 



C The list of M^') functions 



Here we give the explicit form of the M^, j) , M^, M^, M 2 ( g ) and functions (see 

(2.26)). 



(i = l;j = l): 

M 3 (11) = -2i<9 (a 7 CT X2i) + 2zdi(ao7rX2i), 



(ii) 



2<9 (ao7aX2i) ~ 29i(a 7rX2i), 



M 13 = -2«7r9i(aoX2i) + 2i7 CT 5 (Q!oX2i), 
= -47 T <9i(a X2i) + 47 CT <9 (a X2i), 
(ii) _ 



M 



1 
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7<T«o7 CT Xn — 7(7^1 (aoX2i) + -^Cao7rXn 2 ^ 



M. 



(ii) 



J 

4^ 



J 
2^A 



(C.l) 
(C.2) 
(C.3) 
(C.4) 

C<9o(a X2i), ( c - 5 ) 



21 



J 



-7 CT ao7(rXii 



1 1 

r 7<7-9i(aoX2i) - -^Cao7rXii + ^^C^o(aoX2i)- (G.6) 
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(i = l;i = 2): 

M 3 (12) = -2^ (ao7aX2 2 ) + 2^(0^X22), (C.7) 

M| 12) = 2d (a o7(T x 22 ) - 2di(a 7rX22), (C.8) 

M^ 2 } = -2i lT di{a X22) + 2*7^(00X22), (C.9) 

M[l 2) = -4 7r a 1 (aoX22) + 47^0(00X22), (C.10) 

M (i2) = _'t/A laaQl(7Xl2 _ ?!^ 7(T ai(aoX22) + -^=C"o 7 rXi2 - ^=C<9o(a X22), (C.ll) 



M ^i) = _ W^^ aQ ^~ i2 _ ^^^^(ao^) - _^^ ao7rX ~ 12 + ^^C^o (0:0X22). (C.12) 
(i = 2;j = 1): 

M 3 (21) = 2i9o(ao7<xXii) - 2zdi(o 7rXii), (C.13) 

Mf l) = -2do(a 7aXii) + 2<9i(o 7tXii), (C.14) 

M x (21) = 2i 7r 9 1 (aoXii) - 2i 7(7 ao(a Xii), (C.15) 

M x (21) = 4 7T ai(a Xn) - 4 7CT a (a Xii), (C.16) 

M (2i) = _1v/^ 7ctQ , o7ctX2i + ?!^ 7(T Q 1 ( aoXll ) + _^^ ao7rX21 + _i_^ ( Q!oXll ) ) (C.17) 

M (2i) = _f|^ 7CTQ , o7CT ~ 2i + ^^ 7(T g 1 ( ao ^ 11 ) _ J_ ( A ao7rX ~ 1 _ _l_^ ( 9 ( Q , oXll ). (C.18) 
(i = 2;j = 2): 

M 3 (22) = 2^0(007^X12) - 2idi(a 07r xi2), (C.19) 

Mf 2) = -2do(a o7(7 Xi2) + 2ai(Q 07r Xi2), (C.20) 

M x (22) = 2i 7r a 1 (a Xi2) - 2i 7(7 a (a Xi2), (C.21) 

M[f ] = 4 7T ai(aoXi2) - 4 7(j ao(aoXi2), (C.22) 

M 23 2) = -^j^7aao7aX22 + ^y^7 CT <9i (00X12) + -^C«o7rX22 + ;^C<9o(ooXi2), (C.23) 

M ^(22) = _11^ 7ctQ , o7ct ~ 22 + ^^ 7(T Q 1 ( ao ^ 12 ) _ _!_^ ao7rX 22 - ^-^C<9o(a Xi2)- (C.24) 
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With the choice of the v^Ui fields (2.36), it is easy to show the following equivalence relations 
between the equations (2.27) - (2.30): 

M| 12 ) + M{f = <^ 



fj 21) + M x ( 3 21) 

4 22) + M{ 3 22) 








„ J21) 

M\ ' 


„ ,(211 


= 0, 


, ,(221 

M\ ' 


, ,(221 


= 0, 


n 11 
M\ ' 


(111 

- M^ 1] 


= 0, 


m] 12) 


C1 9"! 

- -^14 


= 0, 






= 0, 




-m| 2 ) 


= 0, 






= 0, 


M§*> 




= 0. 



Mg X) - = <=! 

Mg 2) -M£ 2) =0 <^ 

m; 3 2i) - m 2 (2i) = o <^ 

Ml (22) - M 2 ( f = <=} 

These relations show that the independent equations are the ones gives in the formulas (2.37) and 
(2.38). 

D Useful identities 

We collect here some useful identities used throughout the text. These identities can be proved 
by direct verification, with the use of the equations of motion in the form (B.l) and (B.2). 
The following identity is equivalent to the constraint (2.31) 

- fart + v> 2 >2 - r 2 ^ 2 = -idi {w) . (d.i) 

The identity (2.40) can be derived by utilizing the following identities: 



fate + r 2 ^i - - mi = -~ {r 2 ^2 - ^ + ^Wi - , (d.2) 

V^Vi - i>2^\ - ^i*V4 + ^ii>2 = -—^ (i>*i>i + ipii>i - i>2^2 - . (d.3) 

The constraint (2.41) is obtained by computing explicitly the left and the right hand sides, and 
showing that both are equal to: 

j [41l&M% + ^2>2>1 " tfltflMl ~ V^lVW + 2^2*^2 + 2V4>2V>2>2 



2?'\/A 5?vA 

+2^*^1^2 + 2^^^] - —j- [V4Vi + ^> 2 ] + ~yt [rM2^iM*2 
+V' 2 >iV'2^iV'iV ; 2] - —jr WJfaip'iipi + ipiipiip 2 ip 2 - ^2*^2^* ^1 - 1P1M2 ^2] 

+ 16 ' (A) J (ff2)2 W*^v4*^Wi - ^i^MWx] + W^i^i^Vi 

+^1*^2*^2 + ^2^1*^2* + ^[^1^2*] • (D-4) 

The constraint (2.42) corresponding to the zero anomalous term can be proved by first writing the 
left hand side as: 

+ 22 = {ll + l\)C x + (2ihh)C 2 , (D.5) 
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where 



Ci = 3i(a o V ; 2)( 2 «o7T'02) + {2aoj T Tp2)di(a i/j 2 ) + do^o^) (2007^2) + (2ao7aV'2)<9o(aoV'2) 
+ 3i(a V ; i)(2a;o7TV'i) + (2aoj T ipi)di(a ipl) + d (a ipl)(2a ^ a ip 1 ) + (2ao7aV ; i)5o(ao'i/ ; i) 
+ d (a ^2)di(a ijj2) + ^(aoV'^^i^o^) + d (a ^i)di(a ^l) + 9 o (a O '0i) 5 i( a oV ! i) J ( D -6) 

and 

C 2 = (2a 7rV ; 2)5i(aoV'i) + (2a 7rV ; 2)5i(«oV , i) + <9o(a ^2)(2ao7<7^i) + do(aoV'2)(2ao7aV ; i) 
+ (2a O 7r'0i)5l(«o^2) + (2a 7rV , i)^i(«o^2) + d (aoipi)(2a jaip 2 ) + ^(ao^iX 2 ^^) 
+ d 1 (a ipl)do(a ip2) + 9o(«oV ; 2)5i(«oV'i) + d (ao^i)di (0-0^2) + ^(ao^doCaoV'i)- ( D - 7 ) 

Then, one simply shows that C\ = 0, and C2 = 0. 



E Matrix elements of 

In this appendix we give the explicit formulas for the matrix elements of Vt(x). For completeness, 
we give the explicit dependence of these functions on the coupling constants g 2 and g 3 (see footnote 
9). 

= 2 J 2 ' X3X4 ~ X4X3 ^ 4 j3 X2X3X4X4> 



"12 = 8J3 1X2X3X4X3 - X1X3X4X4J ) 

"13 = -%+ 2J2 (X2X3 - X4Xlj + 8J3 (X2X3X4X2 + X1X2X4X4) 

i\/Ac/2 / , . , \ 
= 2 J 2 ' X2X4 X3Xl ^ 

+ 8J3 (2X2X3X4X1 - 2X1X2X3X4 - X1X3X4X2 + X1X2X4X3) , 

i\f\g 2 ( , ,n n/Xcfe / 

"22 = 2J2 (X3X4 - X4X3J ^JsT X1X3X4X3) 



= 2J 2 (-X1X3 - X4X 2 J 
8J 3 



+ Z g A 73 g3 (-2X1X3X4X2 + 2X1X2X4X3 + X2X3X4Xi - X1X2X3X4) , 



n . . W A 92 ( , . 1 \ . W A g 3 , , 1 \ 

52 24 = + (-X1X4 + X3X2J + 8J3 (~XiX3X4Xi - X1X2X3X3J > 



i\/\g 2 , , n , »vA , 

= 2J2~ ' XlX2 ~ X2Xl ^ 473 ^^2X4X2) 

J2 3 4 = 8J3 (X1X2X4X1 - X1X2X3X2J ) 



E.6 
E.7 
E.8 
E.9 



S2 44 = 2J2 (X2X1 - X1X2J ^3— X1X2X3X1- (E.10 

We note that Q(x) is a symmetric matrix, 

0^=%. (E.ll) 
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Moreover, its elements satisfy the following properties under complex conjugation: 

J7 X1 = —^33, f2 12 = ~ ^34, ^13 = —^13) ^14 = ~^23, ^22 = — ^44, ^22 = ~ ^44- (E.12) 

F The Dirac structure 

In this appendix we present the Dirac structure of the AAF model. For completeness, here we 
also give the explicit dependence on the coupling constants g<i and 53 (see footnote 9). 



^11 = 1^92 (-X1X2 + X2X1) + ^53XiX2X4X 2 + (F- 1 ) 
+ ^Ji 9 2 (X1X2X2X3 + X1X4X1X2) + 

+ JJ&92 (-X1X2X3X1X2X4 - 3x1X2X4X1X2X3) > 



(0 ^22 = 1^92 (-X1X2 + X2Xi) - ^53X1X2X3X1+ (F.2) 

+ ^Ja92 (-X1X2X1X4 - X2X3X1X2) + 
A? 

+ JJ&92 (3xiX2X3XiX2X4 + X1X2X4X1X2X3) > 

(^ _1 )l2 = 53 (-X1X2X3X2 + XlX2X4Xl) + (F.3) 

+ ~^Jl9\ (-X1X3X1X2 + X1X2X2X4 - X1X2X1X3 + X2X4X1X2) - 

A A^ 

- ^55253X1X2X3X4X1X2 + ^Jq92 (X1X2X3X1X2X3 - X1X2X4X1X2X4) > 

(O- 1 )" = "^52 (X1X3 + X4X 2 ) + (F.4) 

+ ^53 (-X1X2X3X4 + X2X3X4X1 + 2X1X2X4X3 ~ 2 XlX3X4X2) + 

+ £745! (-X1X3X2X3 + X1X4X1X3 + X1X4X2X4 - X2X4X2X3) - 

- 4J652 (X1X2X3X1X3X4 + X2X3X4X1X2X4) 1 



(^" 1 )23 = ^52 (xsx'i + X2X 4 ) + (F.5) 



+ ^53 (X1X2X4X3 - X1X3X4X2 ~ 2X1X2X3X4 + 2 X2X3X4Xl) + 

+ ~^Ji92 (-X1X3X1X4 + X2X3X1X3 + X2X3X2X4 - X2X4X1X4) + 



.,2 

l\1 



+ JJq92 (X1X2X4X2X3X4 + X1X3X4X1X2X3) 
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(jri )l3 = i + £^ ff2 ( X2X ' 3 _ X4X ;) + (f.6) 

+ £j^53X2X4 (X1X4 - X3X 2 ) + 

+ £^#2 (X1X2X3X4 - X1X3X1X3 + X1X3X2X4 - X1X4X2X3 - 

- X2X3XlX4 - X2X4XlX 3 - X2X4X2X4 + X3X4XlX 2 ) + 

+ £^5253X1X2X3X4 (-xix' 3 + 2 X 2 xl) + 

.,3 

+ 776^2 (-X1X2X4X1X3X4 + X2X3X4xix 2 X 3 + 2x1X2X3X2X3X4 - 2x1X3X4X1X2X4) - 



4J 6 
3iA 2 g 2 4 
4J 8 

V 

2J 2 



XlX2X3X4XlX2X3X4> 



(0- 1 ) 24 = i - l —^g 2 (X1X4 " X3X 2 ) + (F.7) 



+ 53X1X3 (X2X3 - X4Xl) + 

+ -^92 (X1X2X3X4 - X1X3X1X3 - X1X3X2X4 - X1X4X2X3 - 

- X2X3XlX4 + X2X4X'lX 3 - X2X4X2X4 + X3X4XlX' 2 ) + 
+ 7J55253X1X2X3X4 (-2X1X3 + X2X4) + 

• x 3 

+ jjed (X1X2X3X2X3X4 - X1X3X4X1X2X4 - 2xiX2X4xix 3 X 4 + 2x2X3X4X1X2X3) - 



4J 6 

3^A 2 <7| 1 1 1 1 

4J8 XlX2X3X4XlX2X 3 X4^ 



^33 = 52 (-X4X3 + X 3 X4) + ^J3 53X4X2X4X3+ (F.8) 



+ £^1 (X1X4X4X3 + X4X3X2X3) + 

+ JJq92 (-X2X4X3X1X4X3 - 3X1X4X3X2X4X3) , 



^44 = ^7^52 (-X4X3 + X 3 X4) - £^53X3X1X4X3+ (F.9) 
+ £^5 2 (-X 2 X 3 X4X3 - X4X3X1X4) + 

.,3 

+ 176^2 (3X2X4X3X1X4X3 + xiX4X 3 X2X4X3) , 
(Jr 1 )^ = ^J753 (-X4X1X4X3 + X3X2X4X3) + (F.10) 
+ 7745! (-X4X3X1X3 + X 2 X 4 X4X3 - xix 3 X4X3 + X4X3X2X4) - 



4J 4; 

£^5253X4X3X2X1X4X3 + ^Jq92 (xix 4 X3XiX4X3 - X2X4X3X2X4X3) , 



3 

iX , , i\z 



33 



We note that the matrix Q is also symmetric: 

(n-%- = (Q- 1 )^. (F.n) 

Furthermore, its elements satisfy the following properties under involution: 

Or 1 )^ = -(fi- 1 )!!, (Q'Tu = -(^ 1 )22, (^- 1 )^ 4 = -(^ 1 )l2, 

(n-^u = -(^~ 1 ) 3 2, = -(^- 1 )i3, (n- 1 )^ = -(^ 1 ) 24 . (F.i2) 

Finally, we address a subtlety related to the validity of the Jacobi identity. Even when consid- 
ering the simplest forms of the Jacobi identity, such as: 

{Xl OMi)) ,{xi [x<r(2)) ,Xl K(3))}} = 0, (F.13) 

where P c stands for all cyclic permutations of (1,2,3), one arrives at meaningless expressions, since 
the left hand side fails to vanish. One has to remember, however, that the correct way to define the 
Poisson brackets in field theory is done by introducing functionals as in (4.27): 

F(x) = J dt;c(x,0xi(0, (F.14) 

where c(x, £) is some smooth generalized function with some properties on the boundary. The 
corresponding Jacobi identity is then satisfied, provided some general conditions on the functionals 
(see [3,45] for more details). 



G Computational details for the algebra of Lax operators 

In this appendix, we collect some useful formulae necessary for the derivation of the algebra of 
Lax operators in section 5. There we also use the following representation for the two-dimensional 
Dirac matrices: 



a 



1 

-1 



1 






1 



(G.l) 



The functions Q (x; /x) and \x; fj) used in the decomposition (5.1) of the spacial component 



of the Lax connection are: 



£o CT) = 77 [-X3Xi + X4X2 - X1X3 + X2X4] 



4J 
1 

4J2 



[X2X3X4Xl + XlX3XiX2 + X1X2X4X3 ~ X1X2X3X4] 



£l = gj ^3X1 + *4X2 + XlX3 + X2X4J 



+ 



" Vj 

il.f 



-■I + j^J (X4Xl - X3X 2 + X1X4 - X2X3) + (-X1X3 + X2X4) 



X2 + 77X2X3X1 + 7j (X2X4X2 - X1X3X2) 



1 



+ 



Xi - 7JX1X2X3 + 7j (x2X4Xi - X1X3X1) 



16J 2 

1 

16J 2 



X1X2X3X4X2 



X1X2X3X4X1 



(G.2) 



(G.3) 



(G.4) 
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16J 2 



1 



16J 2 



1 



XlX2X3X4X3 3 



XlX2X3X4X4 



(G.5) 



where we dropped the dependence on x and on the spectral parameter [i to avoid cluttering. 
The functions fn, and appearing in the expressions (5.14)-(5.15) and (5.16) are: 



Til = Xl {x'3^13 ~ X4^23 1 + X'1^33 1 ~ x'2^44) + X2 ("X^ll + X4 fi 24 ~ x'l^ 3 4 + xk^u) 

+ X3 (Xa^ll - X4^12 + Xl^M - xfoll) + X4 (-X3^2 + X2^22 ~ Xl^23 + X^m) > ( G - 6 ) 
= 2^ [x2X3X4X3 fi ll + (-X1X3X4X3 ~ X2X3X4X 4 ) ^12 

+ (X2X3X4X1 - X1X3X4X2 - X1X2X4X3 + X1X2X3X4) ^13 + XlX3X4X4^22 1 

+ (-X2X3X4X1 + X1X3X4X2 + X1X2X4X3 - X1X2X3X1) fi 24 ~ XlX2X4Xl^3 3 1 

+ (XlX2X3Xl + X1X2X4X2) ^34' - XlX2X3X2^44 1 ] ~ 4 X2X3X3X4^11 

+ 4 (X1X3X3X4 + X2X4X3X4) fi 12 + 4 (X1X4X2X3 - X2X3XlX4) ^13 

+ 4 (X2X3X1X3 - X1X3X2X3 + X2X3X2X4 - X2X4X2X3) ^14 ~ 4 XlX4X / 3X4^22 1 

+ 4 (xiX3XlX4 - XlX4XlX3 + X2X4XlX4 - X1X4X2X4) ^23 + 4 (X1X4X2X3 - X2X3X1X4) ^24 

~ 4XlX4XlX / 2^33 1 + 4 (XlX3XlX2 + X2X4XlX2) ^34* ~ ^X2X3XlX2^lh ( G - 7 ) 

T 12 = 4 X2X3X4X3^U - 4 (X1X3X4X3 + X2X3X4X4) ^12 

+ 2 (2X2X3X4X1 - X1X3X4X2 - 2 XiX2X4X3 + X1X2X3X4) ^13 + 2 (X1X2X3X3 - X2X3X4X2) ^ll 
+ 4 XlX3X4X4^22 1 + 2 (-X1X3X4X1 + X1X2X4X4) ^23 

+ 2 (-X2X3X4X1 + 2X1X3X4X2 + X1X2X4X3 - 2X1X2X3X4) ^24 - 4 XlX2X4X'l^33 1 

+ 4 (XlX2X3Xl + X1X2X4X2) ^ - 4 XlX2X3X / 2^44 1 - ( G - 8 ) 

The expression that appears in (5.21), after taking into account the explicit expressions for the 
Dirac structure, has the form: 



In this appendix, we list the functions N- ^2) appearing in the matrices Ni(x,y; (ii, (jq), 

i = 0,1,2, which provide the non-ultralocal decomposition of the Lax algebra for the fermionic 



2r n + 




2i (x3Xi + XiX2 + X1X3 + X2X4) 




H The List of (/ii, /i 2 ) functions 
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Wadati model (5.26). 





M2) 




M2) 


<Vi 


M2) 




M2) 




M2) 




M2) 




M2) 


Jvf( Ml 


M2) 




M2) 




M2) 




M2) 



i2(M2j L**4(A*l)X3 + Z 3(Ml)X 4 J + 73 X 3 ^ X4> 



7 AJ 
2 

2 

2 

2 

2 

/X7 



J 2 



, / s r ., / s / , / n /1 , 2a 2 (/i 2 ,^i) // Ufcipv, Hi) /, 
i2(M2j L-"4(/*l)X3 - *3(Ml)X4j + 73 X 3 ^ X4> 



jl 



J 3 



M^l) F4(/^2)X3 + <3(^2jX4j H 3 X3 3 X4> 



J2 



jl 



, f u., / s / x 1 , 3ai(/i 2 ,^i) , 3^1(^2,^1) / 

£2(^2) F4(Ml)X3 + t 3 (^l)X4j H 3 X3 3 X4> 



Jl 



J§ 



t2(M2j h«4(Ml)X3 - <3(/^l)X4j + 3 X 3 3 X4> 



J2 



, / w / s , i s 1 . "2(^1,^2) , ^2(Atl,M2) , 
*2 (Ml) [-«4(At2)X3 - £3(M2jX4j H 3 X3 3 X4> 



jl 



Jl 



1 

7! 
1 



[-Q!l(/ii,yU 2 )X3 + i/?l(/il,/i 2 )X4] , 

[ai(/i 2 , iii)x 3 - i/3i(M2, Ml)X4] , 



J 3 
1 

7l" 



r [«2(M2, Ml)X3 - ^2(^2, Ml)X4] , 



-a 2 (Mi> M2)X3 + ifhifii, M2)X4] , 
-^3(^1)^3(^2) - i£*4(A£i)<*4(A£2)] ■ 



(H.l 
(H.2 
(H.3 
(H.4 
(H.5 
(H.6 
(H.7 
(H.8 
(H.9 
(H.IO 
(H.ll 



Here ai(/Ui,/i 2 ) and /3j(/ii, /u 2 ), i = 1,2 are the following functions of the spectral parameters: 



ai(/ii,/i 2 ) = -/ 2 (/ii)/ 3 (^ 2 ) + [2 - <i(Mi)] £4(^2) 
a 2 (/xi,M2) = h{^i)h{^2) + [2 + £*i(/ii)] ^4(^2), 
/8i(A*i,M2) = '2(^1)^4(^2) + [2 + ii(A£i)] '3(^2), 
feifJ-i, ^2) = -h((J>i)h(fJ>2) + [2 - £*i(mi)] £"3(^2) 



(H.12 
(H.13 
(H.14 
(H.15 
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